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ABSTRACT 


The major contributing factor towards the development of modem 
aeroengines is use of a high gas entry temperature. Thus, an intensive cooling of 
turbine blades is essential. Film cooling using discrete jets in conjunction with the 
internal convective cooling are often deployed to keep the blade temperature under 
the metallurgical limit. The art of blade cooling has established itself as a very 
important subject in aerospace industries. The successful implementation of blade 
cooling guarantees the long life of an aeroengine. The determination of optimal 
configuration of the turbine blade cooling requires extreme care which calls for 
highly accurate computational methods. The flow field near the coolant holes is very 
complex due to the interaction of coolant jets with the hot cross flow. The aero- 
thermal behaviour of the flow depends on a large number of parameters such as 
coolant emission geometry, blowing ratio, density ratio and free stream conditions. 

The present study aims to illustrate the variation of complex aerodynamics 
and heat transfer phenomena which occurs in a jet-cross flow interaction with the 
changes in injection angle, blowing ratio and jet exit velocity profile. The time 
dependent mass averaged three-dimensional compressible Navier-Stokes equations 
are solved based on finite volume formulation. The explicit second order accurate, 
cell centered, time-split algorithm of MacCormack is used. Even on the modem 
computational platforms, the numerical simulation of film cooling applied to turbine 
blades needs a considerable CPU time. Thus, the present study is limited to the film 
cooling on a fiat plate, which meaningfully illustrates the physics of jet-cross flow 
interaction. Admittedly, effects of rotation and surface curvature are absent in this 
analysis. The turbulence is simulated by the algebraic Baldwin-Lomax model. It is 
worth examining the performance of the Baldwin-Lomax model for this particular 
simulation, considering both simplicity and economy of computation. The Baldwin- 
Lomax model is modified with a relaxation technique in an attempt to account for 
the upstream turbulence history. It is worth noting that the Baldwin-Lomax model is 
capable to capture the complex flow physics reasonably well. The inclusion of 
relaxation technique improves the accuracy of the predictive procedure. 
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Chapter 1 


INTRODUCTION 

1.1 GENERAL 

Quite a few components of the gas turbines used in modem aircrafts, operate 
in extremely harsh environment where the ambient temperature routinely exceeds the 
melting point of available alloys. The materials that are usually available cannot stand 
temperatures in excess of 1300 K. The need to protect/cool the solid surfaces exposed 
to high temperature is an important task. At the same time, the maximum temperature 
of the cycle can be increased by cooling such surfaces. Incorporation of various 
cooling techniques and development of new materials have brought about tremendous 
increase in turbine inlet temperature. Today, the commercial aircraft engines operate 
with the maximum cycle temperature in excess of 1600K. 

Various cooling methods concerning the turbine components have emerged 
during last twenty five years. The most useful method is to introduce a secondary 
fluid into the boundary layer on the surface to be cooled/protected. The various means 
of introducing this secondary fluid into the boundary layer include ablation, 
transpiration and film cooling. In ablation cooling, an added coating or heat shield 
decomposes due to sublimation and other highly endothermic processes. As a result a 
significant quantity of gas enters into the boundary layer. In transpiration cooling, the 
surface is usually made of porous material and the secondary fluid enters the boimdary 
layer through this permeable surface. 



Both ablation and transpiration cooling are primarily designed to protect the 
region where the secondaiy^ fluid enters the boiindary layer. They are highly effective 
since a considerable portion of the heat flux towards the wall can be taken up by the 
injected coolant where the heat load is the highest. In addition, the gas entering the 
boundary layer effectively thickens it, and the heat transfer rate is decreased. These 
two methods, however, do suffer from some disadvantages which preclude their use in 
many applications. The ablating material is not in general renewable and so ablation 
cooling has been restricted to systems with high heat fluxes of short duration. 
Examples are re-entry vehicles. This restriction does not apply to transpiration 
cooling since a coolant can be continuously introduced through the porous surface. 
However, the porous materials usually do not possess high mechanical strength which 
is essential for certain applications like rotor blades of the turbines. Also, small pore 
size often leads to clogging and a resulting maldistribution of the coolant flow. 

In film cooling, albeit a secondary fluid is added to the boundary layer, there 
exist considerable differences in operation and even in goals as compared to ablation 
and transpiration cooling. The fundamental difference is that the film cooling is not 
primarily intended to protect the surface just at the location of coolant addition. 
Rather it protects the region downstream of the injection location. Film cooling is 
thus the introduction of a secondary fluid (coolant or injected fluid) at one or more 
descrete locations on a surface exposed to a high temperature environment. The 
surface covered with a thin film of cool air is not only in the immediate region of 
injection but also in the downstream region. The interaction of the coolant jets with 
the hot mainstream is illustrated in Fig. 1 . 








The geometry and the flow field at the point of injection are significant 
variables in film cooling. The secondary fluid can enter through a continuous slot at 
some angle to the wall surface and the mainstream (Fig. 1(a)). 

In gas turbines the air is usually drawn from the exit of the compressor. This 
involves some penalty in turbine efficiency. However, a considerable gain in turbine 
performance outweighs this loss. In the case of a rotor blade, pressure losses can be 
minimised by expanding the air through adjustable slots/nozzles oriented tangentially. 

Although the injection geometry can influence the film cooling performance in 
two-dimensional study of flow, the effect is usually of second order as compared to 
the geometrical effects in three-dimensional film cooling. In this latter, the injection of 
secondary fluid is not uniform across the span. The cool air is injected at isolated 
locations, often through descrete holes on the surface of intest. This can lead to a 
situation equivalent to jets of secondary fluid being blown off the surface. For 
structural reasons it is usually impossible to have a truly continuous injection slot, and 
therefore, interrupted slots and even rows of multiple holes have been used (Fig. 1(b)). 

Although film cooling has primarily been used to reduce the convective heat 
transfer rate firom a hot gas stream to an exposed wall, it can also be used to shield a 
surface from thermal radiation if the radiation absorptivity of the injectant is high. 
This can be effectively accomplished with gas particle suspensions or a liquid coolant. 

The introduction of secondary fluid into the boundary layer during the film 
cooling may be considered to produce an insulating layer (film) between the wall to 
be protected and the gas stream flowing over it. Alternatively the injected fluid can be 
considered as a heat sink that effectively brings down the mean temperature in the 
boimdary layer. The introduction of the secondary fluid into the boundary layer at a 






temperature lower than that of the mainstream and its resultant mixing with the fluid 
in the boundar>' layer reduce the temperature in the downstream. Injection of the 
secondary gas can increase the boundary layer thickness and the mass of the fluid 
entrained into the boundary layer from the free stream. The increased boundary layer 
thickness tends to decrease the heat transfer to the wall (Figl'2). 

The focus of the present problem is primarily on the applications, like cooling 
of turbine blades at the high pressure section of the turbine. It is needless to mention 
that the surface to be cooled in such applications is a curved one. In addition, the 
rotation wall have an intrinsic influence on the flow field. However, for the simplicity 
of analysis, cooling of a flat surface has been considered in the present investigation. 
Eventhough the situation under consideration reveals an ideal case, the physics of jet 
cross flow interaction and the mechanism of coolant layer formation are well 
understood from such a study. 

1.2 SOME BASIC DEFINITIONS 


Let us consider the flow described by Fig. 1-2 The external stream static 
temperature is T^, the film temperature is T^, and the local wall temperature is T,^,. The 
heat flux between the mainstream and the surface in presence of film cooling is given 
by 

q = h(T„-TJ 

where, T^^^, is called the adiabatic waill temperature and represents the surface 
temperature of a perfectly insulated surface and h is the film heat transfer coefficient. 
In absence of film cooling, the recovery temperature is given by 


T=T 




2 



i+r 




14 ^-^Mi 


In the above expression, r is known as recovery factor and expressed as 

T -T 

f — r ^00 

~ Ull2C„ 

00 p 

It has been found by Eckert and Drake (1972) that r = VP, for laminar flows 
and r = for turbulent flows. For incompressible flows, T^, and To^ are 
nearly identical and the film cooling effectiveness may be expressed as 



where, is the stagnation temperature of the coolant near the entry. For high speed 
flows (M„ > 0.3), the cooling effectiveness is given by 


7 = 



The cooling effectiveness is basically a function of the blowing ratio 
density ratio (Pa/Peo)* the geometry of the injection hole, the 
distance from the injection hole (X), and the injection angle (a). 

The experimental results of film cooling have been well documented by 
Goldstein (1971). Jabbari and Goldstein (1978) have developed a correlation for the 
film cooling effectiveness from their experimental study as 


7 = 


1.9(Pr) 


2/3 


1 + 0.329 


jpoo 


\Cp^j 
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iP2U2lpJ^JL 


\ 

Hi. 

\pj 


Re, 


1 - 0.25 



/9=l + lixRe"^Re. 




sin a 


where W is the molecular weight, Rej = and L is the equivalent slot width 

(characteristic length of the injecting hole). 

We may mention that the early developments of computational methods 
addressing the film cooling problem were aimed at matching the above mentioned 
correlation. However, a tremendous development has taken place in the areas of 
numerical modeling and the experimental techniques ever since the paper of Jabbari 
and Goldstein (1978) was published. A concise description all these studies have 
been documented in a subsequent chapter (Chapter 2). 




Fig 1.3 A Film cooled Turbine Blade section ' 


\ ' 
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1.3 SCOPE OF THE PPJESENT INVESTIGATION 


The importance and application of film cooling have already been enumerated 
in the earlier section. The pertinent parameters which govern the cooling 
effectiveness significantly are the blowing ratio, the angle of injection of the cooling 
fluid, and the velocity profile of the coolant at the exit of the injection hole. A need 
has been felt to analyse and imderstand the influence of the pertinent parameters on 
the effectiveness for a situation with multiple injection holes on a plate. This analysis 
will form a building block which can be enhanced to mimic a wide range of operating 
conditions encountered in practice. The aim of the present investigation is to study 
primarily the effect of injection velocity profile, angle of injection and blowing ratio 
on the film cooling effectiveness. The flow regime for such applications is usually 
turbulent. The analysis concerning nuances of turbulence is beyond the scope of this 
work. However, we have used Baldwin-Lomax (1978) model for the present 
computation. For the purpose of further improvement of the predictive procedure, we 
have embarked on a modification of the Baldwin-Lomax model. Finally an effort has 
been undertaken to compare the numerical results vidth their experimental counterpart 
due to Schiffer and Hennecke (1991). 



Chapter 2 


LITERATURE REVIEW 

2.1 FILM COLLING OF FLAT PLATE GEOMETRY 

A brief review of the literature pertaining to film cooling over a flat plate and a 
turbine blade is presented here to illustrate the state-of-the-art. The research in the 
field started with study of slot injection and porous injection in which adiabatic 
effectiveness over the plate could be obtained by using a two dimensional model of 
the flow field. The study on the flow field due to injection through a single hole and 
a row of holes followed these investigations, The evaluation of complex flow fields 
necessitated a three dimensional analysis. Subsequent investigations, in this field 
include the influence of different parameters such as the blo'wing ratio, density ratio, 
length to diameter ratio of inlet pipe, the free stream and jet turbulence, the state of the 
incoming boundary layer and so on. An excellent review by Goldstein (1971) 
describes some of these factors and earlier experimental investigations on film 
cooling. 

The fundamental research on film cooling started since 1950. But the first 
research report was published from the University of Miimesota by Goldstein et al. 
(1968). They considered injection through a single circular hole for their study. This 
was followed by a report of Goldstein et al. (1970) which addresed the issues 
concerning the injection through a series of holes. 



Till 1970, the study was carried out with density ratios around 1.0 and very 
low blowing ratios. In one of the earlier computational works (Launder et al. (1978)) 
a three dimensional analysis of a turbulent jet from a circular hole in a cross flow was 
accomplished. They considered a low blowing ratio. 

Foster et al. (1980), have conducted experimental study on the flow and film 
cooling effectiveness for the case of injection through a row of holes. The 
experimental results include the influence of various injection angles, at different hole 
spacings and blowing rates. They reported that at small injection angles the 
effectiveness is high for low blowing ratios. On the other hand at high injection 
angles the effectiveness is high for high blowing ratios. Closer hole spacings 
improved the lateral coverage and alleviate jet lift-off effects. 

Demuren et al. (1986) reported a systematic study of film cooling with a three 
dimensional calculation procedure. They presented the effect of different parameters 
on the cooling effectiveness for one row of holes. They used standard k-e model for 
evaluating turbulent stresses and heat fluxes. They deployed the locally - elliptic 
procedure of Rodi and Srivatsa (1980), which allows the full elliptic treatment to a 
small region containing reverse flow. However, in the near and the far fields, the 
model reverts to the partially parabolic and fixlly-parabolic schemes respectively. 

Although Foster and Lampard (1980), Firth and Jones (1988) have covered a 
wide range of density and blowing ratios, one feature they had in common was the use 
of large length-to-diameter ratio of the entry port. Such a configuration is not 
representative of a gas turbine. Sinha et al. (1991) have conducted experiments using 
length-to-diameter ratio of 1.75. They have also calculated adiabatic film cooling 
effectiveness for a large number of density and blowing ratios for a row of 35° 



inclined holes with pitch-to-diameter ratio of 3. Schiffer and Hermecke (1991) 
experimentally studied the influence of gas temperature on the film cooling 
effectiveness. The temperature of the hot gas was varied between 420 K and 780 K. 
The free stream Mach number was 0.4. The experiments were performed on a flat 
plate with a row of inclined holes at an angle of 30° and pitch-to-diameter ratio of 2. 
As per the author's knowledge, this is perhapes investigation, cited in open literature 
where the test facility used hot gas to estimate film cooling effectiveness over a fiat 
plate. Thus the flow environment in the above mentioned experiment closely 
represents the mechanism of film cooling in a gas turbine. 

Amer et al. (1992) presented a comparison of different two-equation 
turbulence models for prediction of film cooling effectiveness for two rows of holes. 
The models used include k-oo model, a modified k-© model, that k-e model together 
with a modification of anisotropic turbulence. The analysis reveals that the prediction 
of effectiveness depends on the blowing ratio and the distance downstream of the 
injection hole. They used the PHONICS code to solve the governing equations. The 
study indicates that the two-equation turbulence models do not work well for the film 
cooling problem especially in the vicinity of the holes. 

Leylek and Zerkle (1994) have conducted a large scale computational analysis 
and compared their results with earlier experimental results of Sinha et al. (1991). 
The three dimensional Naveir Stoke equations have been solved to analyse the jet- 
crossflow interaction. Their analysis compares favourably with the experimentally 
observed, velocity profiles (due to Pietrzyle et al., 1989) at the exit port due to 
interaction with main stream for different blowing ratios. 
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In the case of low length to diameter (L/D) ratio and for blowing ratio of 
approximately 0.5, the exit velocity profile is almost similar to a fully developed 
turbulent profile (l/7th low profile). For blowing ratios above 0.5, the velocity profile 
is skewed with the peak shifted towards the upstream direction. Similarly for blowing 
ratios below 0.5, the peak is shifted towards the downstream. Near the hole, the 
effectiveness is not well predicted. Moreover, the lateral spreading of the coolant is 
underpredicted due to an isotropic nature of the flow field near the injection region. 

The aerodynamics of a coolant jet in an environment, which is representative 
of what prevails in a gas turbine has been numerically investigated for film cooling 
applications by Sarkar and Bose (1995). The time dependent mass averaged NS 
equations coupled with the compressible form of a two equation low Reynolds 
number k-e model have been solved using an explicit finite volume formulation. The 
results present the effect of blowing rates, slot angle, and the free stream turbulence at 
the inlet.* They have also analysed the performance of five different turbulence 
models in predicting the flow-field and heat transfer phenomena created by the jet- 
crossflow interactions pertaining to film cooling. 

Sarkar and Bose (1997) have presented another numerical study of film 
cooling with three dimensional calculation for injection through descrete holes on a 
flat plate. The results presented illustrate the variation of complex aerodynamics and 
heat transfer phenomena in a jet cross flow interaction. The simulation comprises 
solving the three dimensional and Navier Stokes Equations (mass averaged) along 
with the algebraic turbulence model of Baldwin and Lomax (1978). 
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2.2 ADVANCED RESEARCH IN FILM COOLING 


Although the present work is restricted to the numerical simulation of discrete 
jet film cooling over a flat plate, some recent work of film cooling over turbine blade 
are sited below: 

Sarkar and Bose (1995) demonstrated the aero-thermal behaviour of coolant 
jet over the VKI turbine blade for different blowing ratios and angles of injection. A 
three-dimensional simulation by Sarkar (1997) also indicates the mechanism of 
coolant layer formation over a highly curved vane and a layer of contra-rotating 
parallel vertices. This study illustrates that the formation of coolant layer on the 
pressure surface is critical and the regions between two adjacent holes are not well 
protected. 

Garg et al. (1997) investigated the effect of velocity and temperature 
distribution at the exit of the hole on film cooling of turbine blades. Garg and Abhari 
(1997) simulated film cooling on a transonic rotating rotor blade and compare their 
results with the experimental observation. On the pressure side, the experimental 
result is under predicted. The turbulence is simulated by the Baldwin-Lomax model. 
The authors prefer to use the Baldwin-Lomax model considering the fact that, it is 
simple and relatively economical compared to two equation models. They also 
indicate that the two equation models are not found suitable for film cooling 
applications. 
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Chapter 3 


STATEMENT OF THE PROBLEM 

3.1 OBJECTIVES OF THE PRESENT WORK 

The main objectives of the present work are to simulate the complex jet cross 
flow interaction in film cooling applied to a flat plate with injection of secondary fluid 
through a series of circular holes and study the effect of different parameters on the 
effectiveness of film cooling. The parameters considered here for study include the 
blowing ratio, the angle of injection and the velocity profile of the secondary fluid at 
the exit of the blowing hole. 

The Reynolds averaged (Farve averaged) three dimensional Navier Stokes 
equations and the energy equation are solved using a finite volume approach. For the 
modeling of the turbulent flow, Baldwin-Lomax algebraic model has been used. The 
Baldwin-Lomax model has been modified in order to invoke the effects of turbulence 
near the injection hole. Explicit MacCormack method is deployed for solving the 
governing equations. 

As mentioned in the literature review. Schiffer and Hennecke (1991) have 
accomplished film cooling experiments with hot gas as the mainstream over a flat 
plate. Their results for the cases of some specific governing input parametershave 
been considered as the basis of comparison of our predictive procedure. Leylek and 
Zerkle (1994) have indicated about various exit velocity profiles at the injection holes 



for different physical situations. These profiles include the fully developed turbulent 
pipe flow and the skewed velocity profiles. We have considered all these possible 
cases in our study including a uniform profile at the exit of the hole. The results 
obtained have been discussed in a subsequent chapter. 

3.2 GOVERNING EQUATIONS 


The Reynolds averaged Navier-Stokes Equations of fluid motion and energy 

equation can be non dimensionalised using the flowing non-dimensional variables. 

. X . y . z 
X y 2 = — 




All the barred quantities in the above equlities are Favre averaged (1965) 
quantities for a turbulent flow. On the basis of the nondimensional procedure 
mentioned above, the Reynolds number and the prandtl number are defined as 


Pr = ^ 
k 

The coefficient of viscosity (p.) is a function of temperature of the gas and is 
calculated by the Sutherland’s equation in the following way: 


/I = 1.456x10- 


r3/2 


7 + 101.4 



The equations of fluid motion and energy after nondimensionalisation can be 


written as: 


1. Continuity Equation 
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2. X-Com « onent of Momentum Equation 
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3 . Y-Com » onent of Momentum Equation 
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(3.3) 


4. Z-Component of Momentum E » uation 
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5. Eneryy Equation 
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In the above equations, Pr is the Prandtl number of the fluid which has ben 
taken as 0.7 and Pr, is the turbulent Prandtl number. A value of 0.9 has been used as 
turbulent Prandtl number. 

The above set of Favre averaged Navier-Stokes equations and the energy 
equation can be written in the following strong conservative form. The prefix is 
omitted at this stage for convenience. However, all the variables are necessarily non- 
dimensional. 
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3.3 TURBULENCE MODEL 
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The quantity, n, in the earlier equations is the eddy viscosity or the turbulent 
viscosity and it is calculated using an algebraic model called Baldwin-Lomax model 
(1978). The Baldwin-Lomax model is described herein. 

3.3.1 Baldwin-Lomax Model 

Baldwin-Lomax (1978) model is a two layer algebraic eddy viscosity model in 
which p., is given by 
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where, 6 is the smallest normal distance from wall at which the inner and the outer 
formulae are equal. 


For the inner region, Prandtl-Van Driest formulation is used which is given by 



(3.7) 


where 


I = Ky\- exp(-/ / ) 


and 

|ivj is the magnitude of w (vorticity) given by 



dU 


A 2 




3V 


dz) \ dz dx.) 


and 


y* — Pw^ry _ V/^w^K> y 

P. P. 

For the outer region (From Closure Formulation) 
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The quantities and are determined from the function 
F(y) = - exp(-/ / A*)\ 


The Klebanoff intermittency factor. 



1+5/%^] 

Ujif is the difference between maximum and minimum total velocity in the 

profile. 

In effect, the distribution of vorticity is used to determine length scaler so that 
the necessity for finding the outer edge of the BL is removed. 

3.3.2 Modification Of Baldwin-Lomax Model (Relaxation Model) 

A relaxation model is generally used to account for the complicated interaction 
of scaling in spatial response to the turbulent structure in a rapidly changing flow field 
such as flow over a rearward or a forward facing step. 

When a secondary coolant is injected into the hot cross flow, there occurs a 
blockage effect with local acceleration of the crossflow over the coolant jet and a 
wake region is evolved with flow reversals. Thus in the present environment, it has 
been decided to incorporate the relaxation technique of Waskiewicz et al in an attempt 
to account for the upstream turbulence history effect. 

The relaxation eddy viscosity is given by 

forx>x, (3.9) 

where, is the calculated equilibrium value of eddy viscosity obtained from the 
Baldwin-Lomax model and Pup^ denotes the value of the eddy viscosity at a location 
(Xp) upstream of the injection. At this location the surface pressure rise begins owing 
to the injected fluid. The quantity 6^ is the boundary layer thickness at Xp. The 
parameter X is the relaxation length scale which represents the lag-response of the 



turbulent stress to the sudden change of the flow gradient. The condition X = a 
implies that the eddy viscosity at all streamwise locations downstream of the 
interaction is frozen to initial upstream values (HapJ. It is assumed that X = 7.5 will 
characterise the present turbulent mixing. 


Chapter 4 

NUMERICAL METHOD 


4.1 FINITE VOLUME FORM OF GOVERNING EQUATIONS 

The imsteady three dimensional mass averaged Navier-Stokes Equations 
described in the previous chapter are marched in time using explicit, second order 
accurate timesplit, cell centred finite volume scheme based on MacCormack [1972] until 
the steady state is reached. 

Rewriting the equation (3.6) 

dU ^ ^ dG 
— + — + — + — = H 
^ Sc ^ & 

Integration of the above equation over a finite volume leads to its integral form as, 
— \UdQ. + j Q.ds = ^HdQ = 0 (As H = 0) 
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where 

Qp denotes the volume of the cell (ij,k) 
ds denotes the surface normal vector 


(4.1) 



Fig. 4.1 





Q is the flux vector of E, F and G 

These systems of equations can be solved in the x, y, z coordinate system for a 
control volume shown in Fig. 4.1 by application of the Explicit MacCormack algorithm. 
4.2 EXPLICIT MACCARMACK METHOD 

The unsteady equations are solved by LI, LJ, and LK operators as explained by 
MacCarmack (1972). In each of the operator the solution proceeds in one of the three 
coordinate systems, x, y and z. Each operator has the predictor and the corrector steps. 
In the predictor step, the fluxes are calculated using the variables at known time le\ el. In 
the corrector step, the predicted value is corrected to get the value of the variable at the 

present time level. In the predictor step the flux Q.dsis evaluated by backw^ard 
differencing whereas in the corrector step it is evaluated by forward differencing. This 
procedure is repeated till steady state solution is reached. 
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Using the above operators, the ' is obtained by operating the LI, LJ, and LK 
operators in the following sequence. 
u;;' = [U(A/ / 2)LKiAt ! 2)LI{At)LK(At / 2)U{At / 2)]U” 


where, 

At = min {A^, 2 Aty, 2At^} 

The surface vector (Fig. 4.1) used above is written in mathematical form as 

“^1234 ~ ^ ^3) 

F, 3 ,r 24 unit vector of the diagonals. 

The volume of the cell ijk (Fig. 4.1) is calculated mathematically as. 
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4.3 COMPUTATIONAL DOMAIN 


4.3.1 Grid Generation 

The geometry consists of a flat plate with desecrate holes of diameter, D = 4mm. 
An algebraically generated non-uniform Cartesian mesh consisting of 85x40x20 grid 
points is used as the computational domain. The pitch to diameter ratio is taken as 2. 
Because of the use of rectangular grids, the edge of the exit hole has to be approximated 
by steps. Due to limited computational facilities, the grid size is restricted to the 
configuration mentioned above. The grid is generated by placing a number of fine 
uniform grid points within the hole region and stretching the grid in all the three 
directions. A fine exponentially stretched mesh spacing is employed in the near hole 
region, and a course uniform grid is used far away from the hole of injection. 

4.3.2 Boundary Conditions 

Boundary conditions are specified at all the six boundaries of the computational 
domain under study. At the inflow boundary, the total pressure and total temperature are 
specified. The static pressure is extrapolated from the interior and hence the other 
variables at the inlet are updated for each time integration. The exit plane is positioned 
far downstream of the injection (100 D) so that a Neumann condition can be applied for 
all variables at this plane. On the solid surface, i.e. on the plate face, no slip conditions 
are imposed alongwith normal temperature gradient as zero to satisfy the adiabatic wall 
condition. The upper boundary of the computational domain is chosen sufficiently far 
from the plate and a ‘slip-wall’ type boundary condition is applied at this plane of 
boundary. The flow is symmetric with respect to x-y plane through both the centre of the 
injection hole and midway between the two adjacent holes. A symmetric condition is 




imposed at these boundaries by prescribing z-component of velocity' vector and setting 
gradients of all other dependent variables in z-direction equal to zero. 

4.3.3 Injection Conditions 

The injection condition within the hole is difficult to specify. As a part of the 
present study, different velocity profiles are used for the injection condition. Leleyk and 
Zerkle (1994), and Pietryzk et al. (1989) have investigated computationally as well as 
experimentally that the injection velocity profile will be same as fully developed 
turbulent pipe flow in the case of a blowing ratio of 0.5. For the other cases the peak of 
the profile will be skewed from the axis of symmetry. Pietryzk et al. (1989) have 
mentioned that a uniform velocity profile will predict almost similar result to that of 1/7 
th power law. In the present work all the above possible velocity profiles are used as 
injection condition and the effect is observed. 



Chapter 5 


RESULTS AND DISCUSSION 

5.1 GENERAL 

In most of the film cooling applications discrete holes are used for introducing 
the coolant. The individual jets emerging form desecrate holes in a row interact with 
the cross flow and thereby create a strong three dimensional flow near the injection 
region. The jet trajectory is bent by the interaction with the main stream and this 
bending produces a pair of counter rotating vortices inside the jet which modify the jet 
cross sectional shape to a kidney like form. The individual jet interacts with the 
neighbouring jets and the main stream. This interaction creates a layer of counter 
rotating parallel vortices which develop into the cooling film downstream of injection 
The flow topology corresponding to the hole injection is shown in Fig. 5.1 . 

For the present work, the experimental investigation of Schiffier and 
Hennecke (1991) is numerically simulated using the compressible Navier-Stokes 
equations with the Baldwin Lomax model. The relaxation techmque is also 
incorporated into the Baldwin Lomax Model to take into the account of the upstream 
effect. The flow configuration and computational domain are shown in Fig. 5.2. 

The boundary conditions have to be prescribed at all the confining boundaries 
of the computational domain which have already been discussed. However the 




Flow configuration for injection through a row of holes 







specification of the flow condition at the exit of the jet is known to be a critical issue. 
A viable approach is to include the whole part into computational domain [Zerkle and 
Leylek,1994]. Experimental investigation indicates that for a low L/D ratio, which is 
closer to the gas turbine film cooling applications, a fully developed profile is evolved 
at a blowing ratio of about 0.5. At a higher blowing ratio, the jet exit profiles are 
skewed with a peak shifted to the upstream side whereas for a low blowing ratio the 
peak is shifted towards the downstream (Fig. 5.3). It has been decided to examine the 
effect of exit velocity profiles on the film cooling, with a low’ L/D ratio. The 
experiment of Schiffer and Hennecke (1991) has been taken up as the basis of 
comparison in the present invetigation. The L/D ratio used is 3.0 and the flow 
environment is similar to what is found in gas turbine applications. The free stream 
Mach number is 0.4, free stream temperature is 658 K, and Re„ based on the hole 
diameter is 2.5x10*. 

All computations have been performed on a DEC Alphastation - 2000. The 
solution is assumed to be converged, when the value of residual have dropped by 
atleast 4 orders of magnitude. The steady state solution achieved in about 10000 
iterations with the Baldwin Lomax models. The present computations are carried out 
with a 85X40X20 grid with pitch to diameter ratio of 2.0. Figure 5.4 shows the grid 
on x-y plane. The grid is stretched near the wall and near the injection region to 
capture steep gradients of flow variables. The first grid point above the wall is so 
chosen that the magnitude of y* remains below 5.0. Because of the rectangular grid 
used, the edge of the exit hole has been approximated by steps. 
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5.2 THE FLOW AND THE TEMPERATRUE FIELD 


The computed velocity field, downstream of a row of holes along the hole 
centreline plane is presented in Fig. 5.5. Here the effect of turbulence is simulated by 
using the Baldwin-Lomax model (with built-in relaxation technique) and an exit 
profile, approximated by a 1/7 power law. Figure 5.5 indicates a strong jet crossflow 
interaction in the near field region. The velocity profile in this region exhibits two 
zones namely (a) the momentum deficit zone indicating the development of a shear 
layer between the inner layer of the injected jet and the crossflow, and (b) the region 
right above the momentum deficit zone, with a slight velocity over shoot caused by 
the accelerating cross flow. Figure 5.6 indicates the velocity field on the x-z plane at 
a horizontal cross-section very near to the surface. This figure shows that, the region 
beneath the jet is ventilated by the mainstream sideways and therefore a low 
effectiveness may occur just at the immediate downstream of the injection. At this 
angle of injection there is no sign of flow separation at the immediate down stream of 
the hole. A more complex flow structure is observed on the cross-stream planes. The 
cross-stream velocity vectors are presented in Fig. 5.7. The counter rotating vortices 
are well captured and the strength of the vortices are prominent upto a X/D of 1 0. The 
strength of this counter rotating vortices is reduced than beyond a distance of X/D 
equal to 1 5. TThis indicates that the jet tends to mix with the mainstream and the three 
dimensionsality is eventually reduced. However along the centre line between any 
two holes, there is some rolling of fluid and this rolling is even predominent at X/D 
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Fig 5.8 Computational centreline temperature contours downstream 

of a row of holes 
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Fig 5.9 Predicted temperature contours on different cross stream 

planes indicting the jet trajectory and lateral diffusion 


equal to 20. This is what culminates due to the interaction of the coolant layer coming 
from two blues. 

The jet trajectory can also be determined by the temperature field. Figure 5.8 
depicts the temperature field on a x-y plane along the centre line of the holes. The 
non'Uniformity in the thermal field is prominent in the near field region with a jet like 
trajectory . The kidney shaped cross-section of the vortices are well captured through 
the cross stream temperature contours (Fig. 5.9) for X/D upto 10.0. At a downstream 
distance of X/D = 20, the three-dimensionality is ceased and the temperature field 
shows more uniformity (Fig. 5.9(d)). 

5.3 ADIABATIC FILM COOLING EFFECTIVENESS 

Figure 5.10 shows the predicted film cooling effectiveness (r|). The results are 
compared with the experiments of Schiffer and Hennecke where the exit velocity has 
been assumed to be uniform. As seen from Fig. 5.10(a), the distribution of film 
cooling effectiveness along the centre line of the holes compares well with the 
experimental data in the range of X/D > 15. However, the effectiveness is 
underpredicted in the near field of the injection hole (X/D <15). This under prediction 
may be explained, in light of cross stream ventilation effect. The experimental results 
do not reveal such behaviour. Figure 5.10(b) illustrates the variation of effectiveness 
along the centreline between two holes. It is evident that in the far downstream (X/D 
> 30), the computational results match the experimental data reasonably well. 
However in the near - field region, effectivenss is once again underpredicted. 
Presumably the lateral diffusion of the jet is underpredicted by the present turbulent 
model. The pointer from the turbuelence literture indicates that the near field region 
is dominated by anisotrpic turbulence. Any model based on isotropic eddy viscosity 



assumption admittedly fails to capture the actual distorsion of strain rates. The above 
stated fact is more clearly appreciated by the distribution of lateral variations of film 
cooling effectiveness at different cross-sectional planes in the near-field region (Fig. 
5.10(c)) . Figure 5.11(a), (b) and (c) depict the distribution of the film cooling 
effectiveness onthe plate, at the similar loctions (viz., along the centre line f the holes, 
along the centre line between the holes and along different cross-stream planes) as 
seen in Fig. 5.10. The velocity profile of the emanating jet in Fig. 5.1 1 corresponds to 
a 1/7 th power law profile. Interestingly, in all the cases, the predicted results 
corroborate the experimental observations of Pietrytz (1989) in the downstream 
region. Here the trend is same as previous case. 

Figures 5.12(a), (b), and (c) show the comparison between the predicted and 
experimental results for a polynomial exit velocity profile. The polynomial velocity 
profile in this case is so chosen that the peak of the profile is shifted towards 
downstream direction. Such a profile corrosponds to the experimental observation 
with high blowing ratio. The comparison makes it evident that the polynomial exit 
velocity profile with a positive shifts fails to simulate a realistic flow situation. 

Figure 5.13(a), (b) and (c) show a similar comparison between the experiment 
and the predictive procedure using a polynomial exit velocity profile with a shift 
towrads the upstream. It may be mentioned tht this particular situation corresponds to 
a low blowing ratio. However, the comparison depicts the inability of the simulation 
in the case of a polynomial exit velocity profile with negative shift (shift towards the 
upstraeam). 

Sarkar and Bose (1995) have observed that the Baldwin Lomax model can 
predict effectiveness quite well in the far downstream of the injection hole. They also 
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found that the prediction did not match well with the experimental results near the 
injection hole. The} indicated that a relaxation model would be appropriate to predict 
the near hole region. We have made use of the same relaxation model for 
improvement in the predictive procedure. Figure 5.14 (a), (b) and (c) show the 
distribution of predicted effectiveness due to the aboe mentioned relaxation model on 
the plate at different locations. The predicted results reveal a highly encouraging 
comparison with their experimental counterpart. Figure 5.14(a) illustrates that the 
predicted effectiveness along the hole centreline is very close to the experimental 
results for X/D greater than 20. However no appreciable improvement is discerned in 
the prediction of hole centreline effectiveness near the injection hole. Figur 5.14(b) 
shows the distriburtion of effectiveness along the line between two successibe holes. 
The values of predicted effectiveness are somewhat improved near the injection hole. 
Figure 5.14(c) shows the lateral variation of the effectiveness at various cross-stream 
locations. It is clear from this result that the lateral variation of the effectiveness is 
close to the experimental results at X/D = 20. Some improvement is also brought 
about at other X/D locations but the extent of improvement is not noteworthy. 

5.4 PARAMATRIC STUDY 

The jet trajectory which primarily depends on the angle of injection and the 
momentum ratio ,has a remarkable influence on the film cooling. A high trajectory 
implies deeper penetration of the jet into the cross flow, pronounced flow mixing and 
hence less effectiveness. On the other hand low trajectory results in an attached 
coolant layer, less mixing and high effectiveness. The influence of the blowing ratio 
and the angle of injection on the film cooling will be discussed in the following 


section. 



5.4.1 Blowing Ratio 

Blowing ratio, as explained earlier is the ratio of mass flow rate of the coolant 
to the mainstream flow. Figure. 5.15 shows the variation of effectiveness for three 
different blowing ratios. The angle of injection has been taken as 30°. It is clear that 
the effectivness along the hole centre line (Fig. 5.15 (a)) and the effectiveness along 
the line between two holes (Fig.5.15 (b))are high for low blowing ratios at the near 

injection region.This trend is gradually altered as the flow travels downstream. The 
reason may be attributed to the high penetration of jets into the crossflow and 
more sideway ventilation of mainstream as blowing ratio is increased. However, 

in the far down stream, the film formation takes place and the effectiveness 

improves. Fig 5.16 shows the temperature contours in the streamwise direction 
for two different blowing ratios. 

5.4.2 Angle of Injection 

Figure 5.17 (a) and (b) show the variation of effectivenes for three different 
angles of injection. The results are obtained for a = 30°, 60° and 90° with the blowing 
ratio of 0.559. It can be discerned from Figure 5.17(a) that the effectiveness is high 
for the low angles of injection near the injection hole. Whereas for the high angles of 
injection, the effectiveness low owing to the side ventilation effect. It is also clear 
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Variation of adiabatic film cooling effectiveness for three 
different blowing ratios at various cross stream planes 











Fig 5.20 Computational centreline temperature contour, downstream 

of a row of holes for two different angles of injection. 



that the effectiveness in the downstream is almost uniform for high angles of 
injection. The effectiveness along the centerline between the holes is high near the 
hole for the high angles of injection and in the downstream the effectiveness ceases to 
be a function of angle of injection (Fig 5.17 (b)). 

For the high angle of injection, the effectiveness drops immediately after the 
injection to a minimum and increase again followed by a gradual drop. This 
behaviour can be explained by the jet lift off mechanism. The jet penetrates into the 
mainstream and the region beneath the jet is ventilated by the mainstream resxilting in 
immediate drop of effectiveness at the immediate downstream of the injection holes. 
After a short distance, the coolant fluid reaches the wall due to turbulent diffusion and 
thus the cooling effectiveness increases. In the case of low angle of injection, the 
effectiveness increases near the injection hole due to the low jet trajectory. Increase in 
angle of injection leads to a higher jet trajectory, but the jets mixes more rapidly with 
the main stream than that for low injection angles. Figure 5.18 illustrates that the 
effectiveness is almost uniform due to a = 90° at X/D = 10. Hence it can be said that 
at the higher injection angles, umform lateral spreading is achieved at a shorter down 
stream distance from the injection hole. 
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Chapter 6 
CONCLUSIONS 

Film cooling is one of the most efficient means of protecting the turbine blades 
from the direct impingement of the hot gas in the gas turbines. The present 
numerical investigation is devoted towards an in-depth understanding of the 
coolant behavior and its interaction with the mainstream flow. The study has 
been confined to cooling of a flat plate. Albeit being an appro.ximate one. 
the study identifies the relative importance and influence of the governing 
parameters involved in film cooling. The effect of turbulence has been invoked 
by an algebraic model. The present formulation is based on the well knowm 
Baldwin- Lomax model of turbulence. The modification has been incorporated 
to bring about the HisioryEf feds of turbulence quantifies. The effect of 
coolant velocity profile at the hole exit, the blowing ratio and the angle of 
injection on the effectiveness of cooling has been studied for a wide range of 
parametric variations. The effectiveness corresponding to the l/7th power 
law profile compare favorably with the available experimental results. 
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